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We rewrite the imaginary-time formalism of finite temperature field theory in a form that all 
graphs used in calculating physical processes do not have any loops. Any production of a particle 
from a heat bath which is itself not thermalized or the decay and absorption of a similar particle 
in the bath is expressed entirely in terms of the sum of particle interaction processes. These are 
themselves very general in meaning. They can be straight forward interactions or the more subtle 
and less well-known purely interference processes that do not have a counterpart in the vacuum. 
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I. INTRODUCTION 

For equilibrium field theory at finite temperature, 
there are two main methods of performing calculations. 
They are the imaginary-time formahsm where one starts 
out in Euclidean space and analytically continues back 
into Minkowski space at the end of the calculation, and 
the real-time formalism where the calculation is done in 
Minkowski space with real time throughout. Because the 
latter has explicit real time dependence, it is therefore 
more suitable for time-dependent problems. However it 
also has the feature of the doubling of the field degrees 
of freedom so that each field acquired a partner and the 
propagators became 2x2 matrices, therefore in the sense 
that the components used in calculations are scalar quan- 
titiesQ at r = versus matrix quantities in the real-time 
formalism the simple analogy and straightforward simi- 
larity to field theoretical calculation in the vacuum are 
absent. Intuitively, the difference between the calcula- 
tion at zero and finite temperature should only be that 
the latter acquires thermal weights in the phase space 
integrations. This clearly is not the case in the real-time 
formalism because one has in addition to deal with ma- 
trix quantities. In this regard, the imaginary-time for- 
malism resembles much more the zero temperature field 
theory. Having said that it must be stated that we are 
well aware of the Braaten-Pisarski resummation where 
perturbation theory at finite temperature must be rear- 
ranged [l]-||] so that it is not just a matter of thermal 
weights between zero and finite temperature. There are 
various applications of this resummation scheme, see for 
example the above references and also [||-@ . These have 
also been recast into the form of kinetic equations for soft 
particles within a heat bath |pT|-p^ . We will choose not 



*Here we are considering everything in the zero tempera- 
ture field theory as scalar quantities, in other words for the 
purpose of discussion here, we make no distinctions between 
scalar, spinor, and vector fields related quantities. These are 
all scalar in the sense they have no doubling of the degrees of 
freedom in the vacuum. 



to complicate matters in this paper and leave resumma- 
tion out, not even including it partially as in p^-[l^, or 
in other words we will consider only particles typically 
at the same scale as the temperature or higher when re- 
summation is formally not necessary. We will address the 
imaginary-time formalism itself and deal with a number 
of issues on its usage. Because of its intrinsic similarity 
with the vacuum theory, it should be possible to bring 
it into the form where the difference is essentially in the 
thermal weights. It will be shown, however, that there is 
a limit to how far this resemblance will go. And we will 
go even further by cutting all internal loops. 

The imaginary-time formalism is a very nice formalism 
in that by following the established calculational rules of 
replacing the thermal discrete imaginary energy sum by 
contour integrals and analytic continuation, the calcula- 
tions can be done very similarly to the familiar vacuum 
field theory ||l^. The needed thermal distributions will 
be there at the end and repeated self-energy insertion 
along a single line will be sorted out even if that line is 
subsequently cut or put on-shell provided the rules are 
followed correctly. Thus the formalism is very compact 
which can hide many physical processes. For this reason, 
although the formalism is quite superior mathematically, 
it can be very unclear when it comes to finding out ex- 
actly what physical processes are involved in a particular 
calculation. The mathematical advantage becomes the 
disadvantage when it comes to the physics. For each 
given calculation involving medium modification of the 
properties of a particle or the production of one, it should 
be possible to write the result in a form that is a sum of 
all the individual contributing processes. These should 
apparently be the phase space integral of each relevant 
process weighed by a product of thermal distributions but 
this is not entirely true as has already been mentioned 
and will be clear later. 

Another related aspect in thermal field theory that 
has not been categorically pointed out and discussed is 
the much richer possibility of interference when a process 
happens inside a heat bath. Although interference graphs 
are well-known to be necessary both at zero and non-zero 
temperature, they have largely kept their T — form 
when being discussed within the finite T context . 
Thus much of the richness was not obvious and remained 
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hidden within the "simple" zero temperature Feynman 
graphs. The authors are only aware of infrequent men- 
tions of this here and there, see ref. [^,|^ which are 
two of the few papers which, as far as the aspects of the 
possibility of interference is concerned, had gone to some 
depths. This may be the case because the existence of 
spectator particles easily masks any forward scattering 
or similar processes and so ensuring that it is impossible 
to tell whether the latter happened or not. When we say 
forward scattering here, we mean it in a generalized sense 
because this can happen to both a fermion or a boson. 
We will discuss this below, pursue further in the direction 
of ^2 23 1 but within the imaginary-time formalism and 
use thermal QED coupled to a massive vector particle 
which is itself not part of the heat bath as an example in 
Sect. fvi. 



II. THERMAL PROCESS = VACUUM PROCESS 
® THERMAL WEIGHTS, 
ALMOST BUT NOT QUITE 

In thermal calculations, it very often concerns the rate 
of production of some particles from the heat bath which 
are themselves not thermalized or how the medium mod- 
ifies the properties of a particle such as its decay or ab- 
sorption inside the hot medium. When this is done within 
the imaginary-time formalism, one can calculate this or- 
der by order (again we are ignoring resummation here) 
and the formalism will yield formulae involving thermal 
distributions and other expressions. This is fine if one is 
interested only in the answer. If on the other hand, one 
would like to know what physical processes contribute to 
the production of a particle or how another get absorbed 
or stimulated to decay inside a medium, one can express 
the formalism in another more physically explicit way 
but equivalent to the original formalism. To this end, 
we now makes the claim that the differential production 
or emission rate of a particle with 4-momentum k per 
unit volume per unit time from a heat bath, but which 
is itself not in the heat bath, is essentially given by the 
discontinuity of the self-energy as 
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for a scalar boson (the generalization to vector bosons is 
straight forward) and 
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= I{k) + J{k) 



(2) 



for a fermion. Each of these can be explicitly expressed in 
terms of a sum of all the physical processes contributing 
to the production of this particle represented by the func- 
tion I{k) and J{k) on the right-hand-side (r.h.s.). The 
I{k) part, because of its resemblance to what one gets by 
cutting rules at zero temperature [ p4|j25| , is known (see 



below). It is essentially the sum of the relevant phase 
space integrals over the allowed kinematical range of the 
probability for each production process weighed by the 
particle distributions. The existence of the J{k) part is, 
however, not well-known at least not to the same extent 
that it will be expressed later on in this paper. It is very 
easy therefore to assume erroneously that there is only 
the contribution from I{k). This is the case for example 
in the papers (2^,^^ We will clarify this in the following 
sections. 

For those processes that do not involve spectator inter- 
ference graph (see Sect, ^for clarification) and therefore 
whose probabilities can be written as squared modulus of 
the corresponding matrix elements, they are represented 
by 



I(.k) - /" d$p {2T:f 5\-k + s,p. 



x\Mv\^ F-p 



(3) 



Here in the energy-momentum conserving delta function 
for each process V, there are the signs Si which depend 
on whether each 4-momentum pi is incoming Si = + or 
outgoing Si = — . F-p is a product of particle distribu- 
tions for each of the participants that has an entry in the 
energy-momentum conserving delta function 



(4) 
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where we used /*^+^ and Z^^-' to denote Bose-Einstein and 
Fermi-Dirac distribution respectively and Si is a sign for 
this purpose. The measure d^-p represents that of the 
phase space integrations of the process P 



d^p = Y[ 
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(5) 



We have used this I{k) above and will use it again be- 
low for generic representation of this kind of sum of the 
contributing processes of the squared modulus type. 

The other generic function J{k) is for the sum of inter- 
ference processes that involve spectator particles in either 
one probability amplitude of the convoluting pair and the 
emission-absorption or vice versa of the same particles in 
the other amplitude of the pair. The function can be 
represented by 

J{k) - h'^'v dip-p,{2TTf5\-k+Y 
p' •' iev 

x{Mp,M*p, + M*p,Mp,) Fp, Tp, . (6) 

Because in these interference graphs, there are emission- 
absorption of particles of the same momentum in a sin- 
gle amplitude and some of these possibilities have their 
origin in thermal self-energy insertions on the external 
lines, which can also be viewed as generalized forward 
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scattering on the external lines, the phase space integra- 
tion measure d^ip, for such a typical process V in general 
becomes 



n 



Pi 



(2n} 



for a subgroup S of all the external lines of the process 
V' in which each member j in the subgroup has rij num- 
ber of thermal self-energy insertions. We find it clearer 
here to say thermal self-energy insertion since this should 
be familiar to the readers, but we will eventually switch 
to the new meaning of the occurrence of any generalized 
forward scattering on the external line j. This will be 
further explained below in Sect. 0. The above phase 
space integrations only take care of those momenta that 
enter into the overall energy-momentum conserving delta 
function. There are other momentum integrations that 
usually being labeled as loop-momenta. These are rep- 
resented by the other measure dip-p, in Eq. (^) above. 
Here it suffices to state that we do not consider loops 
as such because each loop can be opened up and inter- 
preted as emission- absorption of particles with the same 
4-momentum. This interpretation will allow us to treat 
these loop-momentum integrals as phase space integrals 
so the measure dip-p, can be expressed in general in a 
form very similar to Eq. (|^). These other phase space 
integrations will therefore acquire thermal distributions 
as well. The product of these is T-p/ in Eq. (j^). This is 
different from F-pi which retains the form in Eq. for 
the external lines. 

For the decay or absorption of a particle with momen- 
tum k and mass M not thermalized in the medium, we 
can similarly write 

i Disc n(fc) 
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for a boson and 
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Ii-k) + J{-k) (8) 



Ii-k) + J{-k) (9) 



exp(-fcO/r) 

for a fcrmion with the same mass. We are treating the 
absorption as stimulated decay in the presence of the 
medium, hence the width here T is the width in the 
medium. The thermal factor in the denominators of 
Eqs. (|l]), (^) and (^) are correct because in the limit 
T — > 0, there is no longer a medium to produce any 
particle so dR/d^k — > 0. The width, however, will re- 
main finite and is given now by the discontinuity of the 
self-energy in the vacuum. 



III. SIMPLE EXAMPLES 

The relation of the discontinuity of the self-energy to 
the sum of phase space integrals over each contributing 



process given in the previous section can be readily shown 
to hold for simple cases. We review the case of the one- 
loop self-energy as a simple example. For the production 
of a massive fermion from other massless fermions and 
bosons in a heat bath, the required loop graph is the one 
of the usual fcrmion self-energy. Within the imaginary- 
time formalism, we write it as 



d^p M{k,p) 

(2^p2 (fc_p)2 



(10) 



Here we have not stated explicitly which theory is being 
considered except by using the usual fermion self-energy 
graph, it has to be one with vector coupling between 
boson and fermion- antifermion. The numerator is sim- 
ply denoted by a function M since its details are not 
required here. After converting the discrete energy sum 
into contour integration and performing the latter JTof , 
one gets 



E(fc) 
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(11) 



The vacuum part has also been included in the above 
expression. Now performing analytic continuation and 
taking the discontinuity of the self-energy with the 
delta function ^(feo + P + — p\) for the production of 
the massive fermion in mind, only the second and the 
third term in Eq. ( |ll|) have the right discontinuity. For 
production, we have to set k° — > —k^ to get 

iDiscE(fc)- - y _g_(2^)5(„fco +p + |fc _ p|) 
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(12) 



This can be rearranged to 
iDiscS(/c) 



exp(fco/T) 



1 



d^p' 



i27r)S(+\p^)-:^{2n)6^+Hp'') 



(27r')4 ^ ^ ' (2^')4 

X (27:)^ S'-^\-k+p + p') 

x/(+nPo)/(-HPo)AA(-fc,p) 
dR 

d^k one-loop 



2kr 



(13) 



3 



the form we expected from Eq. (||) . The M is essentially 
just the squared of the vector-fermion-antifermion cou- 
pling, which is all there is for the squared modulus of the 
probability amplitude at leading order. 

One could easily use Eq. ( pi] ) to obtain the discontinu- 
ity for the decay of a massive fermion into a fcrmion 
and a boson in the heat bath. The energy conserv- 
ing delta function to get from Eq. ( [Ill ) in this case is 
5(kQ — V — |fc — p\). Only the first and last term from 
Eq. (^l|) contribute. It is simple to follow the same steps 
and arrive at the form given in Eq. (^) or 

iDisci;(/e) 
~exp(-fco/T) + l 

x(27r)'*(5(4)(fc-p-p') 

x(1 + /(+'(po))(1-/(-Hpo)) m,P) 
= 2MT\ , . (14) 

I one-loop ^ ^ 

Note that this could equally have been written with (5^^) 
with / as given in the previous section instead of (J^^-* 
with 1 ± / . Both of the above explicit examples have not 
the J{k) function part because this will come at higher 
orders. More precisely these contributions come in at 
two-loop and higher orders. 



each has an entry in the overall 4-momentum conserving 
delta function. In other words, we are dividing the self- 
energy in two by cutting through line 1 to m. There 
are of course more than one way to group lines into the 
form in Fig. ^ and therefore different cuts possible on 
that graph. Also at any given order, the self-energy will 
be a sum of graphs of the form of Fig. |^. The full results 
must therefore be a sum both over different graphs and 
cuts on those graphs. This sum will eventually become 
our sum over different contributing processes. 



1 




FIG. 1. A self-energy diagram at any order in the coupling 
can always be rearranged in this form with momentum pi flow- 
ing through line i except the last line m where the momentum 
Pm ~ k + Pi. The two blobs may be complex vertices 

or (m -I- l)-point Green's functions and they need not be the 
same, hence the different shadings. There are of course more 
than one possibility of such an arrangement for each graph. 



IV. GENERAL PROOF 



A. The Case of All Bosonic External Lines 



In this section in order to put up a framework so that 
later sections can be more readily understood, we show 
in our own way that Eqs. (|l|), (H), (H) and (^) hold in gen- 
eral. We saw in Sec. |l| that both the I{k) and J{k) func- 
tions possess the same structure of I{k) at least for the 
external lines included in the delta function conserving 
the overall 4-momentum. That is as far as the part as- 
sociated with lines of the contributing graphs that carry 
part of the total energy- momentum, the two functions are 
the same. J{k) has some extra phase space integrals and 
distributions but these have little to do with the over- 
all momentum conservation and they can be viewed as 
substructures. This will be discussed later on. First the 
main structure of Eqs. (0), (||), (||) and (^, which is the 
structure of the function I{k), will be shown first. 

For a self-energy at arbitrary order in the coupling of 
the theory, it can always be rearranged into the form 
shown in Fig. |l| where momentum pi flows through line 
1, p2 flows through line 2 ... etc and the last line m 
carries the boson or fermion momentum k of the self- 
energy under consideration and the sum of the other lines 
above it so that Pm = k + Pi- ^^'^ blobs 

shaded differently need not be the same in general and 
they can be any graphs from very simple to very complex 
connecting the lines 1 to m. Our aim is to put all these 
lines on-shell so that they become the external lines and 



Our own approach to arrive at the results in Eqs. (|^), 
(^, (^ and (^ is by using the imaginary-time formal- 
ism and performing the contour integrals of the loop mo- 
menta pi,p2 ■ ■ -Pm-i, without loss of generality, in that 
order. Furthermore, lines 1 to m will all be taken to be 
bosonic and massless for simplicity. The case that some 
lines are bosonic and some fermionic is a generalization 
which does not affect the proof provided equations to be 
used below are suita bly m odified. This will be touched 
upon briefly in Sect. IV B. Since the discontinuity of the 



self-energy contains a number of different processes, we 
will aim at getting only one process, say that with the 
overall energy-momentum delta function 



S<■^Hk + J2^^P^) 



(15) 



where the s^'s is a fixed set of signs associated with this 
particular process depending on which momentum pi is 
outgoing or incoming. We adopt the convention that 
s = -f for incoming or absorption of a particle and s — — 
for outgoing or emission. This applies also to k. By con- 
centrating on Eq. (p^), any pole in the contour integra- 
tions that does not contribute to this particular chosen 
process need not be picked up and will be dropped from 
the discussion. For each line, which one of the two poles 
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can be picked is determined by the signs Si 's so that there 
is only one pole per line that will contribute to Eq. (15). 
Starting from the Pi contour integration, it can either 
pick the pole on line 1 or that on line m. In either case, 
the P2 integration will have also two possibilities, that 
of line 2 or line m in the former case and line 2 or line 
1 in the latter. Since each contour integration will have 
two possibilities to pick a pole so that there will be 2™~^ 
terms that contribute to Eq. (|l^) after all m—1 contour 
integrations have been performed. Picking a pole on a 
line in a contour integral is equivalent to putting that line 
on-shell and giving it a factor of a thermal distribution. 
Since there are m—1 integrations and m lines, all but 
one line will remain off-shell at the end. This line will be 
cut and put on-shell as well when the discontinuity of the 
self-energy is finally taken. Therefore although in general 
the numerators of the 2™"-'^ terms will be different after 
the integrations because different sets of external lines 
may have been put on-shell, they will become identical 
once the discontinuity has been taken. We take the nu- 
merator A/'(fc,pi,p2, • ■ • ,Pm-i) as the part that does not 
include the thermal distributions and is evaluated at the 
mass-shells of these particles. Because it is the same for 
every term, it can be taken out of the following discus- 
sion as a common factor. We can concentrate entirely on 
manipulating the products of thermal distributions into 
the desired form. Those readers not interested in the de- 
tails of this manipulation could accept on faith that the 
final form of the product of distributions is correct and 
skip to section B to only examine the identities and the 
simple examples in section C. 

After the discontinuity has been taken, there will be 
2™"^ terms, each with a different product of thermal 
distributions. Writing the propagators in the convenient 
form 



1 



for i = 1 to m 
1 



1 
2^ 



E 



POi - SiPi 



1 and 



2 

POn 



P7, 



2p„ 



E 



(16) 



(17) 



where pi = \pi\ for i = 1, . . . ,to — 1 and pm = |fc -|- 
^™~^ Pi\. If the Pi contour integration picks up the pole 
on line 1, there will be the factor —(1/2 + /'^-'(sipi)) 
which includes also the vacuum part. On the other 
hand if it is the pole on line m that is picked, then 
there will instead be (1/2 + /'^•'(— s„iPm)). In the first 
case, subsequent integration of P2 will yield either the 
-(l/2-f/(+)(s2P2)) or(l/2-H/(+)(-siPi-s™p™)) ther- 
mal factor depending on whether the pole on line 2 
or line m is picked. In the other case, there will be 
-(l/2 + /(+)(s2P2)) or (l/2-h/(+)(-sipi- Wm)) from 
line 2 and line 1 respectively. Continuing like so, there 
will be different products of thermal distributions. Start- 
ing from the case where each contour integration picks 



the pole of its own line, then only one of the lines 1 to 
TO — 1 picks the pole of line m instead of its own line, 
next two lines do not pick their own lines and so on until 
none of the p^ contour integration pick their own pole at 

P° 



SiPi. Using the simplifying notations, 

/W=/(+)(,,p,), 



(18) 



(19) 
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I rii-\-i-\-j-\ 



/im-j-i — ■ = f-^'i-SniPm " S^P^ " SjPj ) , (21) 



= /^ + ^ { + SmP,n + S^p, + SjPj +•••), (20) 



- f( + ) 



and 



r^-i/2+/r' 



(22) 



etc. The sum of thermal factors from the 2™"^ terms 
with all common factors taken out can therefore be writ- 
ten as 



m — 1 m — 1 m—1 



m— 2 m—1 771— 1 



(+) 

-m — > , . I 



-E E ( n r)/^!-..-.,/ 



(+) 



m— 3 m — 2 m—1 m—1 



E E E ( n nf^-ra^ 

j=l k=j+ll=k+l i^j^k^l 
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X f 

■'-m-1-2 (m-2) • 
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(23) 



The key to finding a way out of this seemingly endless 
sum of terms is the observation that each term must have 
a partner. By that we mean for every term containing 
fm-i part of the thermal weight, there must be an- 
other one that differs from this term only in this fac- 
tor by having —f^2,_-^_2 (m-2) place instead. 

This is true because no matter which sets of poles were 
picked in the contour integrations, the last integration 
of must be able to pick either its own pole on line 

TO — 1 or the pole on the other line that now carries the 
energy This could be the line to if it has not 

yet been touched or another line with the energy 
shifted there via an already performed contour integra- 
tion. Therefore every term must have either f^^li or 
— f^^"^ -I „ , n\ by virtue of the last contour inte- 

J— m— 1 — 2 (m — 2) 

gration p^-i- They can thus all be paired. Now using 
the delta function Eq. (Ill), we get the identity 
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Jm-1 J-m-1-2 (m-2) 



A + ) _ A+) 
Jm-1 J -m-1-2 



-k"/T 



(m-2) 

(+) fi+) 

m+l+2H |-(ni-2) 



l) fm-lf, 



(24) 



After this is applied to every pair, the common factor on 
the r.h.s. of Eq. (||) can be taken out of Eq. (||). We 
now have 2" 



terms to sort out. 



Shifting focus now onto the thermal weight "-"^ 
— f -I „ , For similar reason as before, each 

■'— m— 1 — 2 (m — J) ' 

term must have either one of the two and these again 
form pairs because of the second last contour integration 
ofPm-2- The identity 
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J m-2 J -m-1-2- 



A+) 

(m-3) ) J |-(ni-2) 



A + ) A + ) 

J m-2J h(m-3) 



(25) 



is a more general form of Eq. ( pij) and can be used to 
remove the last factor in the now common thermal factor 
on the r.h.s. of Eq. (|2^) to get the new common factor 



,-fc"/T _ 1^ A+) A+) A+) 

^) Jm-l Jm-2 Jm+1+2^ |-(m-3) 



(26) 



The number of terms has now been further reduced to 
2™"^. Iterating this thermal factor reduction process, 
all 2™~^ terms can be grouped together eventually into 
one common thermal factor 



X • 



- 1^ f(+^ f(+^ f(+^ 

f J m J m—1 •> 7n — 2 J m — 3 

A+) A+) A+) 

' J3 J2 Jl 



(27) 



The first factor will be divided out by the denominator in 
Eq. (||) with fc" — > — fc" for production or that in Eq. (||) 
for decay. The remaining thermal distributions will be 
f'^^Hp) for absorption with s = + in the delta function 
in Eq. (|l5|) or —(1 + /^^•'(p)) for emission when s = 
— as expected. Combining the sign factors siS2 ■ ■ ■ Sm 
originating from the propagators in Eq. (18) but have 
been left out of the discussion so far, this is the F-p in 
Eqs. d) and (H) 



S1S2 ■■■Sm 



(28) 



Now we turn to the phase space integrals in Eq. (Q). 
From the m — 1 loop integrations, there are already the 
3-niomentum measures (Ppij {2itY for i = 1 to m — 1. 
This together with the 1/ {2pi) from the propagators give 
(i$ in Eq. (^ from i = 1 up to the m — 1 entry. The last 
entry m can be gotten by introducing the identity 



/|^(2.)3^(3)(fc + £.,p,).l 



(29) 



Multiplying this by the remaining 1 /2pm and the energy 
delta function (27r)(5(fc° + ^™ pf) from the discontinuity 
gives 



i2n)S^^-\pl)6('\k + Y,s.p,) 



(30) 



By grouping everything else not discussed above which is 
mainly the numerator N{k,pi,p2, ■ ■ ■ ,Pm) together, this 
will be LMpp if in this particular process the two blobs 
in Fig. m contain no internal loop momentum that runs 
completely inside the blob. If either one or both contain 
internal loops then this will be the product of (A^pA^p + 
A^pA^-p) and T-p in the function J{k). The latter T-p 
factor will carry essentially the hidden loop momentum 
integrations and there will be associated thermal weights. 
We will however adopt an approach in which there will 
not be any loop. But this and the remaining parts of the 
function J(fc) are the subjects of Sect. ^ 



B. The Case of a Mixture of Bosonic and Fermionic 
External Lines 

For the more general case where there is a mixture of 
fermionic and bosonic lines amongst the lines 1 to m in 
Fig. 1^, the proof is somewhat more complicated because 
bosonic distribution can turn into fermionic distribution 
and vice versa when they are evaluated at a pole depend- 
ing on whether this has an imaginary part with a total of 
an integer n or a half- integer n -I- 1/2 of (27rz) imaginary 
energy. For a bosonic self-energy, there must be an even 
number of intermediate fermion lines and for a fermion 
self-energy, there must be an odd number of such lines. 
In the latter case in order not to have a contradiction 
in our momentum arrangement in Fig. ^ the last line 
m must be fermionic whereas in the bosonic self-energy 
there is no such constraint. We could nevertheless ar- 
range all fermion lines to be those at the bottom and 
leaving all boson lines at the top in either case. That is 
if there are n boson lines n < m, we arrange the lines 
so that line 1 to n will be bosonic and line n -|- 1 to m 
will be fermionic. Then for a fermion self-energy the last 
contour integration will give the pair combination 

Jm-l ^ J-m-1-2 (m-2) 

Jm-1 ' J-m-1-2 (m-2) 

^ ^{'^ + 1) /m-1 /,lj+l+2+--- + (m-2) (^1) 

if the line m — 1 is fermionic or 

_A+) ,A-) 

Jm-l ^ J-m-1-2 (m-2) 

^ _A + ) _ A-) 

Jm-1 •/-m-1-2 (m-2) 

^ ~{'^ + 1) /rrJ^l /rn+l+2+--- + (m-2) (^^2) 

if this line is bosonic. This is always true because the 
number of fermion lines connecting the two blobs in Fig. m 
is odd. Here we have used the notation /j*" '^ = 1/2 — . 
For a bosonic self-energy, there is the combination 
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f(-) _ f(-) 

Jrn-l J-m-1-2 (m-2) 

, /•(-) 

Jm-1 I •'-m-1-2 (m-2) 

= ~ -*-) frn-l /m+l+2+--- + (m-2) ' ('^'^) 

The thermal factors can be sorted into the desire form in 
Eq. using the identities 

-{fr-f^)fi:ht'fr^ (34) 

{fr + f':f)fi;hft'fr. (35) 

-{ft^ + f'^fl^hft'f^^' (36) 

iteratively as similarly done before in the previous sec- 
tion. The phase space integration measure and the rest 
arc the same as discussed there. The proof of the main 
form of the function I{k) or J{k) is therefore complete. 
In Sect. 1^, the integrand and substructures in J(fc) and 
J{k) or equivalently the internal structures of the blobs 
in Fig. |l| will be the subject. 

C. Examples 

In this subsection, we give two examples of the thermal 
factor reduction from 2™~^ terms down to one discussed 
in the previous section. First we consider a bosonic self- 
energy with m = 3 and one-fermion loop. In accordance 
to our discussions above, we push all fermion lines to the 
bottom so line 1 will be boson and line 2 and 3 will be 
fermion. Without stating explicitly what process to get 
from the self-energy, we choose the very general energy 
conserving delta function 6{ko + sipi + S2P2 + S3P3). Af- 
ter the discontinuity has been taken, the sum of thermal 
distributions is 

+/r'/i^L3-/^7i^L3). (37) 

This can be reduced by pairing as follows. 

F'= ™(/r^-/i^')(-r'+&3) 

= .1.2.3(/}^^+/i-30(/^'-/il^3) 

= -siS2S,fi+^ft^ft\e-'^"/^-l) (38) 

As another example, we consider a fermion self-energy 
with m = 4. This time only line m is fermionic. The 
sum of thermal distributions now is 

p' 1^ f (+) f (+) f (+) I f (-) f (+) f (+) 

r -.1.2.3.4(^-/1 J2 J 3 +7-4/2 h 

j_ fi-) fi+) , f (+) fi+) fi-) 

~rJl J-I-4J3 ' Jl -h J -1-2-4 

_fi-)fi-) f(+) _ f{-)f(+)A-) 

J -4 J-I-4J3 ■1-4 J2 J -1-2-4 

_fi+)fi-) f(-) I f(-)f(-) f(-) ^ 

Jl J -I-4J -1-2-4 J -4 J -I-4J -1-2-4 J ■ 

(39) 



Pairing as before gives 

F' = .1.2.3.4 (/?' - fV) {fT' - fV-4) 

x{fr-fV-2-4) 
= . 1.2.3.4 (/l+^+/i-4^)(/i+^+/i7U) 

-{ft'+.f'-t2-4) 
= .1.2.3.4/1+^ - 1) . (40) 



V. THE INTERNAL STRUCTURES OF 
I(K) AND J(K) 

The previous sections showed how the discontinuity of 
the self-energy could be arranged into phase space inte- 
grations over all the external lines weighed by their re- 
spective thermal distributions. There we called any line 
that showed itself explicitly in Fig. |l| external and these 
all have entries in the overall energy-momentum conserv- 
ing delta function in Eq. (^ and Eq. (H) . In this section, 
we will deal with the remaining structures in these equa- 
tions. These structures come essentially from the blobs 
in Fig. |l| which have not been discussed yet. 

In the simplest case, each blob consists of a tree graph 
with all lines leading to the external lines 1 to m and the 
line with 4-momentum k. This can be just a few lines 
at low orders or a very large tree at high orders. All 
internal lines' energies-momenta are fixed completely by 
the external lines and there is no additional integration 
and thermal distribution other than those of the exter- 
nal lines. If the two blobs are identical, then there is 
automatically the jA^pp. If they are not, then there is 
the interference M-pM^ + M^M-p with jTp = 1. In 
the case that M^, and M-p are of the same order in the 
coupling, there must also be the possibilities lA^pP and 
|A4-pp so that a larger squared modulus can be formed 
\M-p + Aipl'^ . For familiar field theories such as 4>^, 4>*, 
QED, QCD etc., tree graphs with the same number of 
external legs are of the same order so this grouping into 
a larger squared modulus is always possible. This case 
for J{k) is rather trivial but has to be stated first for the 
sake of completeness before we move on to more com- 
plex cases. So non-identical blobs with no internal loop 
are the simplest examples of interference graphs in J{k) 
that can actually be grouped into the class of function 
J(fc). It is the more complex interference graphs gen- 
uinely belonged to the J{k) function part that are our 
main concern. Towards this we now turn. 

In the previous paragraph, we discussed the case of 
the tree graphs. This type of contributions can generally 
be interpreted straight forwardly as scattering processes 
or other interactions. For cases of a few particles in the 
initial state, there will always be the vacuum counter- 
part. These are quite familiar physically. We need not 
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say anymore on these. What we find less obvious and 
less well-known is the physical interpretation of the case 
that the blobs Fig. |^ have internal loops. Within the 
imaginary-time formalism, it is not too difficult, albeit 
tedious, to perform the calculation of these graphs and 
it is all too easy to lose the feel of what are the actual 
physical processes involved in any particular calculation. 
The nature of these graphs made it even harder because 
of the fact that first there is no vacuum counterpart and 
second they are the non-intuitive interference type. If 
one is not clear in one's thinking and clings on to the 
over-familiar vacuum picture, they could even appear to 
be impossible. 




FIG. 2. An example graph for a blob with one loop. 

Let us consider the simple case that one blob has one 
internal loop and the other none, the procedure in arriv- 
ing at the form of I{k) has already put all the external 
lines on-shell. So the calculation of this blob is equiv- 
alent to the evaluation of a one-loop n-point function 
where n < (m + 1). The rules of the imaginary-time 
formalism stipulated that the energy contour integration 
of the loop must pick out every pole in the propagators 
in turn of all the particles along the lines that form the 
loop. This step is equivalent to putting in turn every 
particle propagating around the loop on-shell. For each 
particle around the loop, there are two poles. One can 
be interpreted as its emission into the heat bath from 
one vertex and its absorption from the bath by another. 
The other is absorption from the heat bath now at the 
first vertex and its emission back into it from the second. 
An example graph for a blob with one internal loop and 
three external lines is shown in Fig. El This graph can 
be turned into the six graphs in Fig.^ after the energy 
contour integration put the three internal lines in turn 
on mass-shell. These are indicated by the now on-shell 
large dotted lines in Fig. ||. 




(ul) 



(u2) 





FIG. 3. Performing the energy contour integration is equiv- 
alent to putting the internal lines in turn on-shell. (ul) and 
(u2), (vl) and (v2), and (11) and (12) are graphs resulting 
from putting the upper, vertical and lower internal line re- 
spectively on-shell. Each pair of graphs corresponds to the 
two poles from each line. 

In our special case, it is easy to turn the blob with the 
loop into a sum of tree graphs with emission- absorption 
of particles of the same momentum. If one now takes a 
step back and looks again at the whole picture, there is 
also the other blob that helps make up the whole self- 
energy which consists of only a tree graph. That for our 
example graph in Fig. |^ will be that shown in Fig. ^. 



FIG. 4. The associated graph to Fig. |2[ 

Remembering that Fig. ^ and |^ together form an inter- 
ference contribution, opening up the loop in Fig. ^ would 
seem to render the interference impossible because there 
must be the same incoming as well as outgoing particles 
in both graphs. However, graphs representing interac- 
tions in a heat bath are not quite the same as those in 
the vacuum. It is a common practice to use the same 
graphs in both situations but it must not be forgotten 
that in a heat bath graphs should be understood to be 
implicitly accompanied by spectator particles that made 
up the multiparticle system. So opening up the loop in 
should be accompanied by a change of Fig. ^ into 
The momentum carried by the spectator line in 
is of course different when this is convoluted with 
each of the three (u), (v) and (1) pairs in Fig. ^. 
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FIG. 5. This is how Fig. ^ should be properly represented 
once the loop in Fig. ^ has been opened up. 

So in this simple example, there will be six terms in the 
function J(fc). Each will have a phase space integration 
of either s = ± form 

dfv'-j0^i2^)S^'Hl'-^') (41) 

originating from the previous loop integration. Because 
these correspond to the emission and absorption of par- 
ticles of the same momentum within a single graph, they 
do not play a part in the overall energy-momentum flow 
and have not an entry in the corresponding constrain- 
ing delta function. The emission and absorption of these 
particles must come with thermal distribution. Because 
of the unusual nature of the emission-absorption or the 
absorption-emission in the same graph, each of these in 
our example with only one internal loop acquires the dis- 
tribution 

Tv'^l/2 + f(+\\f\) (42) 

for boson or 

Tv> = l/2-f^-\\f\) (43) 

for fermion in this simple case. For more complex cases 
where the blobs are not identical and each has several 
internal loops, after opening up each loop in the blobs 
and adding the corresponding spectator on the other blob 
the phase space integration becomes 

d^V'=t[-^,{^^)S^'^\ll-ml), (44) 

i 

where L is the combined total number of internal loops 
between the two blobs. If there are any self-energy in- 
sertions hidden in a blob, then this loop will have to be 
opened up too. Therefore each self-energy insertion will 
be turned into the emission and absorption of a particle 
of the same 4-momentum from a single line. If the line is 
fermionic, this is the usual forward scattering, but there 
is also the case that a bosonic line emits and absorbs or 
vice versa a fermion of the same 4-momentum. Thus we 
label these generically as generalized forward scattering. 
In this case, the internal loop turned phase space integra- 
tions will be of the form of Eq. (0) . The factor of thermal 
distribution when there are many loops becomes 

J-p, =n(l/2 + sj(^-)(|/°|)). (45) 

i 

Because there must be four or a larger even number of 
graphs obtainable from each loop, the sum in Eq. (^ over 
V' is a sum over the main processes as well as over the 
possible subprocesses of emission-absorption of particles 
in each main process. 



The steps in arriving at the form of the J{k) part con- 
tribution to the production from or the decay of a particle 
in a heat bath are now shown. So expressing a produc- 
tion or a stimulated decay of a particle in a heat bath in 
terms of all the contributing processes should not be re- 
stricted only to processes expressible in terms of squared 
modulus amplitudes. In other words, there is the very 
important J(fc) interference contributions in addition to 
the better known I{k) part. In the next section we will 
illustrate all these in an explicit example. But before we 
do that it must be mentioned the vacuum parts still have 
to be regularized and renormalized in the usual way. We 
have assumed that this was implicitly understood. 

VI. A MORE COMPLEX AND COMPLETE 
EXAMPLE 

In this section we will use QED coupled to a massive 
vector particle as an example. The relevant lagrangian is 

C = -^F^'F^, - eA^ - gV^) i, (46) 

where F^" and G^^ are the field tensors for the and 
vector fields respectively. The coupling g is taken 
to be much weaker than e, e g, so only higher order 
corrections in e will be considered. Our heat bath will 
consist only of leptons and photons at a temperature T, 
and the vector is not itself thermalized. A massive 
vector with a mass M ^ T will be sent into the heat 
bath to determine the medium modification of its width 
up to two loops. 




FIG. 6. The leading contribution to the decay of . 

At leading order (LO), there is only the decay into a 
dilepton pair contribution to Eq. (M) coming from the 
one-loop self-energy shown in Fig. This is of course 
only a contribution of the type of I{k). From Eq. (jl^) 
after some substitutions, this contribution is 

x(l-/(-)(po))(l-/(-)(p;,)) lA^l^o- (47) 

Here the numerator is just the familiar Dirac trace of the 
lepton loop, one gets after summing over final spins and 
averaged over initial spin 
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I-^Ilo 



3 



P-P + 



M2 



(48) 



At the next-to-leading order, there are only three self- 
energy graphs but there are many physical processes hid- 
den within these graphs. In the following graphs, the 
thick wavy line is for the massive vector while the thin 
wavy line is for photon. The lines overlaid by the vertical 
dashed line are the exposed external lines of Fig. |l|. First 
we examine those arrangements of these graphs with no 
internal loops within the blobs. By putting the three 
internal lines in the two-loop graph with an internal self- 
energy in Fig. 1^ on-shell, three physical processes emerge. 
We get a sum of the squared modulus amplitudes of two 
compton scattering in Fig. ^, of one decay with photon 
emission in Fig. and of one vector-photon fusion in 
Fig. These are clearly all contributions to the I{k) 
function part in Eq. (^). There will of course be the 
other contributions where the internal self-energy is on 
the lower lepton line. These can easily be taken care of 
by a simple factor of two by symmetry. 



FIG. 7. Two-loop self-energy with internal self-energy in- 
sertion. AH three internal lines in the middle will be put 
on-shell as external lines therefore the two blobs are tree 
graphs. 




FIG. 8. Contributions to compton scattering from Fig. 0. 

2 




FIG. 9. Decay with photon emission from Fig. |^. 

2 




FIG. 10. Vector-photon fusion contribution from Fig. |^. 

The other two-loop graph is the one with a vertex 
correction drawn in Fig. ^ Putting all three internal 
lines on-shell, three physical, albeit interference, pro- 
cesses emerge. They are again compton scattering in 
Fig. ^ three-body decay in Fig. |l^ and fusion in Fig. |l^. 
Because they are all interference contributions, it appears 



that they belong to J{k). As we said in Sect. ^ blobs 
that have no internal loop or have only a tree structure 
could be regrouped with other contributions to form a 
larger amplitude. Examining the pairs in Fig. |^ and [l^ . 
Fig. 1^ and ^ and Fig. |l^ and |l^ and also not forgetting 
that the two other lepton lines together with the pho- 
ton line in Fig. |ri| could also be put on-shell as well to 
give similar contributions, this can of course be done as 
is well-known. We merely broke down the contributions 
into parts so that they could be clearly seen, within the 
framework that we are presenting the paper, which con- 
tribution came from which diagram. Loosely speaking 
the contributions from Fig. |ri| should be in J{k) but it 
is better to tighten the definition so that this function is 
restricted to the genuine, not well-known in-medium in- 
terference contributions to the modification of the width. 




FIG. 11. Two-loop self-energy with vertex correction. 
Putting the three intermediate lines on-shell and there will 
not be any internal loop. 




FIG. 12. Interference contribution to compton scattering 
from Fig. |l^. This can be regrouped with Fig. ^ to form a 
larger squared modulus of a single amplitude. It does not 
genuinely belong to the J{k) class. 





FIG. 13. Interference contribution to decay with photon 
emission from Fig. This can also be regrouped with Fig. ^ 
to form a larger squared modulus of a single amplitude. This 
is not really a J(fc) class contribution. 



FIG. 14. Interference vector-photon fusion contribution 
from Fig. |ll] Again together with Fig. ^ this can form a 
large squared modulus of a single amplitude. Therefore it is 
in the I(k) class contribution. 

After combining the contributions from Fig. Q and |l|, 
we get for the three processes 



2Af r| 

(2^ 



compton 



(27r)(5W(p2 



(27r)(5(+)(p'2) 
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= 2 



X (2^^)^^^^ (2^)''^^'^ {k+p-p'-q) 

x/(-Hpo)(1-/(-'(Po))(1 + /^+'M) 



on 7 



(49) 



decay 



(2^)4 
(2^ 



(27r)4' 

(27r)(5(+)(g2) (27r)45('*)(fc-p-p'-g) 



X (1 + /(+n9o)) (1 - bo)) (1 - K)) 



:|^ldecay 



(50) 



2MT 



fusion 



(g.(2.).W(,.)|V(2.).(-)(,'^) 
X ^ (27r)5(+) (g2) (2^)4^(4) (k + q-p-p') 

x/W(<Zo)(l-/(-'(Po))(l-/'-'(Po)) 
x|A^lLio„- (51) 



emission and absorption of particles with the same 4- 
momentum within one graph. One must of course put in 
the associated spectator particles in the graph originated 
from the blob on the left-hand-side (LHS) of Fig. |^ in 
order for the new graphs to make physical sense. These 
are now depicted in Fig. |l^ and In Fig. 16, the mas- 
sive vector decays into a real and virtual dilepton pair. 
The virtual one then either absorbed a photon from and 
emitted it back into the heat bath or emitted a photon 
into the bath before absorbing one back. In Fig. |l^, the 
virtual lepton either annihilates with one from the heat 
bath to recreate another dilepton pair or it is put on-shell 
via a photon exchange with a lepton in the bath. The 
lepton line that has just been put on-shell is shown in 
dashed line in the figure. While all these are happening 
on one graph, there is merely the vector decay in the ac- 
companying graph. Here there is a photon or a lepton 
spectator. 




FIG. 15. In this two-loop self-energy with internal 
self-energy insertion, only two lines will be put on-shell. One 
of the blobs will have an internal loop with two lines. 



The averaged over initial spins, summed over final spins 
matrix elements for the three processes are 



4e 
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I Icompton 



^.,2 ^ 8e^ff^ fp-p' + k-p' ^ p-p' + k-p 

r 1 decay ^ 
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2ip-pr 

ip ■ q){p' ■ q) 



\M\ 



2 

fusion 
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(52) 



(53) 



(54) 



These amplitudes were calculated in |^ for studying the 
change in Z boson properties in the quark-gluon plasma. 
Note that if one is only interested in the medium modifi- 
cation to the width, it is necessary to subtract the leading 
order and next-to-leading order vacuum contribution in 
Eqs. (0) and 

Now we turn to the examples of one of the main sub- 
jects of this paper. In Fig. |^ instead of putting the 
three intermediate lines on-shell, one could do this to 
two lepton lines without touching the photon as shown 
in Fig. |l^. This results in the blob on the r.h.s. in Fig. |l| 
to have an internal loop made up of two lines. As we dis- 
cussed already in Sect. they will be put on the mass- 
shell in turn by the energy contour integration. Thus it 
is possible to unfold the diagram into several ones with 




FIG. 16. Two genuine J{k) class contributions coming from 
Fig. |l^. Photon with the same 4-momentum is emitted and 
absorbed within a single graph. The same photon is required 
as a spectator in the other graph. 



+ 




FIG. 17. Opening up the fermion line in Fig. |l^ results 
in the absorption from and emission back of a lepton into 
the heat bath or vice versa. These are different examples of 
generalized forward scattering on the fermion line. The lepton 
line from the self-energy loop is now shown in dashed lines to 
distinguish them from the other leptons. 

These interference graphs can be worked out to be 

2Mr| 



= 2 



(Fig. 16) 

^(-l)(2vr)<5'W(p^)-^(2.)<5(+)(p-) 



(27r) 
d'^q 



■(27r)(5(+)(q2) {27r)^S^'^\k-p-p') 



(27r)4' 

x{l-f^-\p,)){l-f^-\p'„)) 

X (1/2 + /W(go)) {M M* + >(*X)(Fig.0) , (55) 
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and 



2Mr 



(Fig. 17) 



x(l-/(-)(po))(l-/(-)(p^)) 

X (1/2 - f^-\lo)) {M M* + >(*X)(Fig. . (56) 

Note that the mass-sheU constraining delta function of p 
in these equations has a derivative because of the double 
pole due to the double propagator (1/p^)^. The sum 
over final spins and averaged over initial spins convoluted 
amplitudes are 

{MM*+M*M)^Fis. 16) 
16e2g2 



[p-p' + {p-q){p' ■ q) 

Vp"^ + 2p ■ q 



1 



p^ — 2p ■ q 



)}, (57) 



and 



{M M* + M*M)i^Fis. 17) 



2ip-l){p-p')-p\p' -l)) 



1 



1 



2p-l+p'^ 2p • ? - p2 



(58) 



Similar to what was done to Fig. one can also put 
only two lepton lines of Fig. ^ on mass-shell as shown in 
Fig. |lj. In this case, the blob on the r.h.s. of Fig. |l| has 
an internal loop consists of three lines. Converting the 
photon to a real one results in the photon absorption- 
emission interference contribution in Fig. p^. I nstead of 
a forward scattering on the lepton in Fig. |16[, now dif- 
ferent fermions participate in the emission and absorp- 
tion. Next either lepton line within the loop can be put 
on-shell. One gets interference contributions where the 
vector particle gets absorbed by a lepton before a dilep- 
ton pair is radiated off at the end and others where the 
vector fuses with a virtual photon radiated from a ther- 
mal lepton to form a dilepton pair. These are shown in 
Fig. All these interfere again with the less interesting 
simple decay graph with a spectator. 




FIG. 18. Only two lines will be external lines in this 
two-loop self-energy and therefore one of the blobs will have 
an internal loop made up of three lines. 




FIG. 19. From the internal loop of Fig. |l^, one also 
gets a genuine interference contribution with photon absorp- 
tion-emission within a single graph. 



— <J 




FIG. 20. A similar real interference contribution to Fig. |l^ 
from Fig. ^ Note that the former internal but now on-shell 
lepton line is different from that in Fig. ^ There is another 
similar contribution from Fig. ^ by opening up the other 
lepton line in the loop. 

The mathematical expressions for the contributions in 
Fig. |l^ and ^ are similar to those in Eqs. ( |55| ) and 

(^) except 1) {—l){2Tr)6'^^\p'^) is replaced by the usual 
(27r)(5'^'*')(p2) and 2) the convoluted amplitudes are of 
course not the same. These work out to be 



(MM* + M*M)^Fis. 19) 
16eV f _ {p-p'f 1 
3 I {p-q){p'-q)\ 



(59) 



and 



{MM* + M*M)iv,^. 20) 
A:e'^g^ p' ■ I ( p ■ I — p ■ p' p ■ I + p ■ p' 



P ■ t r 
p-l I 



+ f4^}- (60) 



3 p - l Ik ■ I — p ■ p' k ■ I + p ■ p' 
The full result up to 0{e^) correction is 



F 



2F, 



LO ^ ^ -L comptoii 
+ 2 r(Fig. 16) + 2 F(Fig.|l^) 



F 



(Fif 



2 F 



(Fif 



(61) 



The factor of two for the various contributions to the 
width is due to the two possibilities of lepton and an- 
tilepton to take part in whatever interactions that give 
rise to that particular contribution. Note that in the 

rdccayi Tfusion 

and Fpig. Iigl, each has linear infrared di- 
vergence enhanced by thelBose-Einstein distribution but 
these cancel in the sum of the three p7| , p8| . 

By opening up the internal loops, each term in Eq. ( |6l| ) 
is now a contribution from a clear physical process and 
not just a vague one-loop and two-loop contribution 
which cannot be readily associated with an interaction. 
In the general case of self-energy graphs at high orders 
so that there are multiple internal loops and overlapping 
lines, the imaginary-time formalism is just too compact 
for the physics to be transparent. In any case if one is 
concerned with the physical processes, a lot of terms in 
the self-energy will not contribute and drop out once the 
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discontinuity has been taken. Therefore there will be a 
lot of fruitless labor, so instead of calculating hoops and 
loops of virtual particles, it is much simpler to disentangle 
all of them and calculate instead the amplitudes of on- 
shell particle interactions. We have shown here that these 
could be both squared modulus processes and purely in- 
terference processes made possible by the presence of the 
heat bath. The latters are not so well-known and it is 
easy to be influenced by the accustomed vacuum picture 
to believe that there are only the former contributions. 

Using these techniques we have, beside working out 
the Z boson width in a quark-gluon plasma also 
calculated the high mass next-to-leading order dilepton 
production from such a QCD plasma. This is reported 
in I29I ] . Some of these have previously been examined in 
for example [T^ . Since we have not included any form of 
resummation, we do not expect the results to agree. But 
there should be an agreement once the resummed version 
of the present work is done. We will leave this as a future 
work. 

As a final remark although we exclusively worked 
within the imaginary-time formalism, the result for any 
of the physical quantities on the LHS in Eqs. (|^), (||), (|^) 
and (|^) obtained within the real-time formalism should, 
of course, be the same. The connection between the 
real- and imaginary-time formalism for the first few N- 
poi nt functions have been worked out by various authors 
[po| , ^ 30 - 3^ . The relevant quantities in our case are the 
imaginary part of the two-point functions. These are re- 
lated to those in the real-time formalism by 

Imn(A:) = -|-tanh(/co/2T) Imnii(fc) (62) 
ImS(fc) = - coth(fco/2r) ImEii(fc) . (63) 

To calculate the imaginary part of the self-energies in 
real-time, one can use any of the m any suggested finite 
temperature cutting rules |^ , ^ ,p3 26 3^,|^,^. Since 
we are more interested in what physical processes ac- 
tually contribute to the production or decay rate of a 
non-thermalized particle than the mathematical rules by 
which one calculates the imaginary part of n(fc) or S(fc), 
we will leave this at that and do not elaborate any fur- 
ther. For the details of the precise connection between 
any aspects of the two formalism, the readers should 
consult the references given here. The message that we 
would like to convey in this paper is that the plwsical 
processes that contribute to the LHS of Eqs. (|l|), (||), (||) 
and (^ are far richer in numbers and stranger than those 
found at zero temperature. The existence of a thermal 
bath permits many purely interference processes to con- 
tribute which have no counterpart in the vacuum. It has 
also been pointed out that it was possible to calculate 
the rates, provided that the finite temperature interfer- 
ence processes were included, by basing on the actual 
physical processes rather than on which N-point or N- 
loop functions. The latter are not immediately physically 
transparent. 
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